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Neutron stars are expected to contain a significant number of hyperons in addition to protons and 
neutrons in tlie liighest density portions of their cores. Following the worlc of Jones, we calculate 
the coefficient of bulk viscosity due to nonleptonic weak interactions involving hyperons in neutron- 
star cores, including new relativistic and superfiuid effects. We evaluate the influence of this new 
bulk viscosity on the gravitational radiation driven instability in the r-modes. We find that the 
instability is completely suppressed in stars with cores cooler than a few times 10^ K, but that stars 
rotating more rapidly than 10 — 30% of maximum are unstable for temperatures around 10^" K. 
Since neutron-star cores are expected to cool to a few times 10^ K within seconds (much shorter than 
the r-mode instability growth time) due to direct Urea processes, we conclude that the gravitational 
radiation instability will be suppressed in young neutron stars before it can significantly change the 
angular momentum of the star. 
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I. INTRODUCTION 

The r-modes (fluid oscillations whose dynamics is dom- 
inated by rotation) of neutron stars have received con- 
siderable attention in the past few years because they 
appear to be subject to the Chandrasekhar-Friedman- 
Schutz gravitational radiation instability in realistic as- 
trophysical conditions (see Ref. ^ for a recent review). 
If the r-modes are unstable, i.e. if the damping timescales 
due to viscous processes in neutron-star matter are 
longer than the gravitational-radiation driving timescale, 
a rapidly rotating neutron star could emit a significant 
fraction of its rotational energy and angular momentum 
as gravitational waves. With appropriate data analysis 
strategies, these waves could be detectable by interferom- 
eters comparable to enhanced LIGO. The r-mode insta- 
bility might also explain the relatively long spin periods 
observed in young pulsars and of older, accreting pulsars 
in low- mass x-ray binaries. 

Recently Jones [||, ^ has pointed out that long- 
neglected processes involving hyperons (massive cousins 
of the nucleons) can lead to an extremely high coefRcient 
of bulk viscosity in the core of a neutron star. Using 
simple scaling arguments he suggests that the viscous 
damping timescales associated with these processes may 
be short enough to suppress the r-mode instability alto- 
gether in realistic astrophysical circumstances. The pur- 
pose of this paper is to investigate this possibility more 
thoroughly. The hyperons only exist in the central core 
of a neutron star where the density is sufficiently high. 
The relevant effects of the r-modes however vanish as r^ 
(where r is the distance from the center of the star). Thus 
the overall effect of hyperon induced dissipation on the r- 
modes depends very sensitively on the size and structure 
of the core of a neutron star. Jones' initial estimates did 
not take properly into account either the structure of the 
r-mode or the detailed properties of the nuclear matter in 
the core of a neutron star. We improve on Jones' analysis 



in several ways: First we evaluate fully relativistic cross 
sections to determine the reaction rates of the relevant 
hyperon interactions. We find that these cross sections 
reduce to the results of Jones |^, ^ in the low-momentum 
limit, but can be about an order of magnitude larger in 
some regimes of neutron-star matter. Second we derive 
new expressions for the bulk viscosity coefficient that are 
appropriate even for a relativistic fluid such as neutron 
star matter. Third we construct detailed neutron star 
models based on an equation of state that includes hy- 
perons and the appropriate interactions among all of the 
particle species present. Due to superfiuid effects the 
temperature and density dependence of hyperon bulk vis- 
cosity turns out to be quite complicated: superfluidity 
increases the viscosity in some cases while reducing it in 
others. And fourth, we use a more accurate model of 
the structure of the r-mode eigenfunction in the cores of 
these stars to evaluate the effects of hyperon dissipation. 

Our analysis shows that hyperon bulk viscosity com- 
pletely suppresses the gravitational radiation instability 
in the r-modes of rotating neutron stars for temperatures 
below a few times 10^ K. We find that the gravitational 
radiation instability acts most strongly at temperatures 
around 10^° K where stars rotating more than 10 — 30% 
of the maximum rotation rate (depending on the details 
of the microphysics) are driven unstable. Our coefficient 
of bulk viscosity is actually several hundred times that 
of Jones 1^, 1^, who suggested that the instability was 
completely suppressed. However, our use of the proper 
r-mode eigenfunction reduces the dissipation by several 
orders of magnitude and we find that there is a window 
of instability. How long it lasts is another matter. If the 
core of the neutron star cools via the standard modified 
Urea process, its temperature remains above a few times 
10^ K for about a day [|. This is enough time for an 
unstable r-mode to grow and radiate away a substantial 
fraction of the star's rotational kinetic energy and an- 
gular momentum into gravitational waves 0. However 



if the core of the star cools too quickly the instability 
might not have enough time to grow before being sup- 
pressed by the hyperon bulk viscosity. The time needed 
for a neutron-star core to cool to a few times 10^ K is re- 
duced to about a second when direct Urea processes are 
able to act [Q, |^. Modern equations of state have large 
enough proton densities in the core that direct Urea cool- 
ing is now expected to act until the neutrons and protons 
condense into a superfluid state, i.e. above about 10^ K. 
The growth time for the gravitational radiation instabil- 
ity in the most rapidly rotating neutron stars is about 
40 s Thus we conclude that the core of a neutron 
star will probably cool too quickly for the r-mode insta- 
bility to grow significantly before being suppressed by the 
hyperon bulk viscosity. 

The organization of the rest of this paper is as fol- 
lows. In Sec. II we provide details of the equation of state 
which we use, including numerical aspects of the evalua- 
tion of various thermodynamic variables and derivatives, 
and the model of hyperon superfluidity that we employ. 
In Sec. II] we present a new derivation of the coefficient 
of bulk viscosity for relativistic neutron-star matter (in- 
cluding several interacting fluids) in terms of the micro- 
scopic reaction rates and thermodynamic derivatives. In 
Sec. [V, we compute the relevant cross sections in order 
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to evaluate the reaction rates for hyperons in a dense 
medium. In Sec. M we combine the thermodynamic ex- 
pressions of Sec. ^IJ with the microscopic reaction rates 
of Sec. [V to obtain expressions for hyperon bulk viscos- 



ity as a function of density and temperature in neutron- 
star matter. In Sec 



FIG. 1: Fermi momenta (in MeV) of the baryons and lep- 
tons as functions of total energy density (in lO^^g/cm'^) for 
Glendinning's equation of state. 



interest to us these equilibrium constraints are 

Up = TT-e +71^ +?^E-, (2.1) 

Aip = - Me, (2.2) 

Mm = Me, (2.3) 

Ms- = Mn + Me, (2.4) 

MA = Mn- (2.5) 

In order to solve these constraints and determine the 
equilibrium state of neutron-star matter, we need explicit 
expressions for the various chemical potentials as tunc- 



V| we evaluate r-mode damping ^j^^^ ^j^^ p^^^^^^ ^^^^^^ densities n, . These functions 



timescales for neutron stars containin g ord inary fluid and 
superfluid hyperons. Finally, in Sec. VII we discuss the 
implications of our results for the r-mode instability in 
real neutron stars, and also attempt to estimate how ro- 
bust these conclusions are. 



II. EQUATION OF STATE 

A. Thermodynamic Equilibrium 

Neutron-star matter is a Fermi liquid which at low 
densities is composed primarily of neutrons n, protons p 
and electrons e. Charge neutrality Up = (where is 
the number density of the i species) and /3-equilibrium 
Mn = Mp + Me (where /ij is the chemical potential of 
the i species) determine the relative concentrations of 
these particles at each density. As the total baryon den- 
sity increases however, it becomes energetically favorable 
for the equilibrium state to include additional particle 
species: first muons /i, and then a sequence of hyper- 
ons S", A, ... These additional particles appear as the 
density exceeds the threshold for the creation of each new 
species. The relative concentrations of the various species 
are determined at each density by imposing charge neu- 
trality and /3-equilibrium. At the highest densities of 



have encoded within them the details of the interactions 
between the various particles in a dense Fermi-liquid en- 
vironment. In this paper we have adopted the expressions 
for these chemical potentials as given by Glendenning's 
relativistic effective mean- field theory 0, [u). Figure |l| 
illustrates the Fermi momenta of the various particle 
species as a function of the total energy density of the 
matter that we obtained with Glendenning's (case 2 [p^) 
expressions for the chemical potentials. Glendenning also 
gives expressions for the total energy density p and to- 
tal pressure p as functions of the particle densities n^. 
These quantities are illustrated and tabulated by Glen- 
denning ]lO| , and we will not reproduce them here. 
Our numerical code reproduces Glendenning's numbers 
quite accurately. 

We are also interested here in some less familiar ther- 
modynamic quantities that are relevant for calculating 
the bulk viscosity in neutron-star matter. These quan- 
tities are easily determined once the full description of 
the equilibrium state is known. In particular the par- 
tial derivatives of the chemical potentials with respect to 
the various particle number densities, a^- = dpi/drij, are 
needed in the expression for the relaxatio n tim e associ- 
ated with bulk visocity as defined in Eqs. ( 3.161 ), ( |3.24 ), 
and ( 3.25| ) below. These aij are easily determined numer- 
ically (or even analytically in some cases) once the full 
equilibrium state is known. Further the thermodynamic 
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FIG. 2: Thermodynamic prefactor 70c — 70 (the difference 
between the "fast" and "slow" adiabatic indices) that appears 
in our expression for the bulk viscosity. 



function, 



7o 



70 



dp dStn 



p drin duB ' 



(2.6) 



appears as a prefactor in the expression for the bulk vis- 
cosity, Eq. ( 3.11 ), that we derive below. Here dp/dun is 
just the partial derivative of the pressure with respect to 
the number density of neutrons (keeping the other num- 
ber densities fixed), and dxn/duB is the derivative of the 
fractional density of neutrons in the equilibrium state, 
in = nn/nB, with respect to the total baryon density ub- 
The left side of Eq. (2.6) has been re-expressed in terms 
of 7oo the "fa st" a nd 70 the "slow" adiabatic indices de- 
fined in Eqs. ( 3.12| ) and ( 3.13| ) below. Figure ^illustrates 
this function for the Glendcnning equation of state. For a 
non-relativistic fluid the pre- factor p(7oo — 7o) is identical 
to a commonly used alternative expression involving the 
sound speeds of the fluid: /3(u^ — Uq) [|2|. However, this 
equality is not satisfied in neutron-star matter. Conse- 
quently itis important to use the correct expression given 
in Eq. (U). 

We have solved the relativistic structure equations for 
the non-rotating stellar models based on this equation of 
state. Figure |^ illustrates the total energy density as a 
function of radius for neutron-star models having a range 
of astrophysically relevant masses. This figure illustrates 
that these stars contain large central cores having mate- 
rial at densities that exceed the and A threshold den- 
sities. The fact that it becomes energetically favorable to 
create hyperons (or even free quarks) above some thresh- 
old density is not really very controversial. However the 
expressions for the chemical potentials as functions of the 
particle densities are not well known, and so the detailed 
properties of nuclear matter at the densities where hyper- 
ons are likely to occur is not well determined at this time. 
This uncertainty translates then to an uncertainty about 
the sizes of the hyperon containing cores of real neutron 
stars. Since the size of this hyperon core determines the 
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FIG. 3: Structure of neutron stars having a range of masses 
using Glendinning's equation of state: total energy density 
(in units of lO^^g/cm^) vs. distance from the center of the 
star. Threshold densities for E~ and A hyperon formation 
are also plotted. These equilibrium structures are computed 
using general relativity. 



strength of the bulk viscosity effects which we evaluate 
here, the implications for the stability of the r-modes are 
correspondingly uncertain as well. 



B. Superfluidity 

Next, we must consider the possibility that the hyper- 
ons in neutron-star matter form Cooper pairs and con- 
dense into a superfluid state at sufficiently low tempera- 
tures. Various calculations are given in the literature of 
the A superfluid gap function Aa ||1^, The A gap 
function is constrained by the experimental data on the 
energ y le vels of double A hypernuclei such as ^^^Be and 
aaB IjlJI, however even so it is probably only known to 
within a factor of two or three. In our numerical analysis 
of the bulk viscosity timescales discussed in Sec. |^ we use 
an analytical fit to the zero-temperature gap function Aa 
as computed by Balberg and Barnea [|l3[. Their calcu- 
lation produces a gap that depends on the total baryon 
density ns of the nuclear matter and on the Fermi mo- 
mentum pf of the A itself. As illustrated in Fig. ^ we 
find that the following empirical fit 



X [0.77 + 0.043(6.2715 



0.88)2] ^ (12.7) 



matches their calculated values for the zero-temperature 
gap energies fairly well. In this expression the gap energy 
is measured in Mev, pp is measured in fm~^, and ub is 
measured in fm~'^. 

The superfluid gaps are not as well determined be- 
cause comparable experimental data on double S~ hy- 
pernuclei do not exist at present. Calculations by Takat- 
suka, et al. |lj] using several models of the nuclear in 
teraction give values of A^- in the range: Aa ^ Aj; 
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FIG. 4: Comparison of the zero-temperature supefluid gap 
function Aa as calculated by Balber g an d Barnea (dots) with 
the empirical analytical fit in Eq. (2.7) (curves). The bot- 
tom and top curves correspond to ns = 0.4 and 0.8 fm~'^ 
respectively. 
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FIG. 6: Temperature dependence of the superfluid gap: ex- 
act values (points) compared to the empirical analytical fit 
(curve) given in Eq. (2.J). The critical temperature Tc is re- 
lated to the zero-temparature gap energy by kTc = 0.57A(0). 
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FIG. 5; Superfluid critical temperatures as functions of the 
total energy density (in units of lO^'^g/cm'^). The solid curve 
is for A hyperons, while the dashed curves are for E~ hyperons 
using either Aj;- = Aa (bottom) or Aj^- = IOAa (top). 



IOAa. We perform two sets of calculations based on the 
extremes of this range. We either set A 2- = A a or 
Aj]- = IOAa. By equality here we mean that the de- 
pendence of As- on and pp (up to the overall factor 
of 10) is given by Eq. (2.7). Using these energy gaps, and 
the Glendenning equation of state, we can evaluate then 
the density dependence of the superfluid gap functions. 
We illustrate these in the form of superfluid critical tem- 
peratures Tc [which are related to the zero-temperature 
gap functions by kT^ = 0.57A(0)] in Fig. |. 

The superfluid gap depends not only on the density of 
the superfluid material, as discussed above, but also on 
its temperature. This temperature dependence will be 
needed to determine the temperature dependence of the 
hyperon bulk viscosity below. The standard BCS model 
calculation fl^ of the temperature dependence of the gap 



is illustrated in Fig. |^. This flgure compares the results of 
the exact calculation with a simple empirical fit to these 
data: 



A(T) = A(0) 



3.4 



0.53 



(2.8) 



Since this fit is quite good, we use it whenever the tem- 
perature dependence of the gap is needed. 



III. BULK VISCOSITY 

Bulk viscosity is the dissipative process in which the 
macroscopic compression (or expansion) of a fluid ele- 
ment is converted to heat. The formalism for calcu- 
lating the bulk viscosity coefficient in terms of the re- 
laxation times of the microscopic processes which ef- 
fect the conversion is well-known — see, e.g.. Landau and 
Lifschitz |]l2j. However, such calculations are generally 
performed for an ordinary fluid like air, in which the 
microscopic processes (typically involving the transfer 
of energy between rotational and vibrational degrees of 
freedom of the molecules) can be taken to be indepen- 
dent. In neutron-star matter there are several relevant 
microscopic processes (involving weak interactions be- 
tween the various particle species), but they are related 
by constraints (such as conservation of baryon number) — 
meaning that we cannot simply use the standard formulas 
for either a single process or multiple processes. Further, 
neutron-star matter is composed of some particles having 
relativistic energies. The standard expressions for bulk 
viscosity are not correct for such materials. In this sec- 
tion we present a modifled and expanded derivation of 
the equations of bulk viscosity appropriate for neutron- 
star matter. 
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Bulk viscosity is due to an instantaneous difference be- 
tween tlie total physical pressure p of a fluid element 
and the thermodynamic pressure p. The thermodynamic 
pressure is determined only by the equation of state for a 
fluid element of given particle number and entropy den- 
sities. It is the value toward which the microscopic pro- 
cesses are driving the physical pressure at any given time. 
The coefficient of bulk viscosity ^ defines the proportion- 
ality of this pressure difference to the macroscopic expan- 
sion of the fluid: 



P-P 



-CV. 



(3.1) 



where v is the velocity of the fluid element. 

Consider now a fluid state that is an infinitesimal per- 
turbation of a time-independent equilibrium state. Let 
Pq and hq denote the pressure and number density (func- 
tions only of position) that describe this equilibrium 
state. To calculate C in terms of the microscopic reac- 
tion rates that drive the system toward equilibrium, we 
re-express both sides of the Eq. ( |3.lD in terms of the 
Lagrangian perturbation of the particle number density 
An = n — uq. Using the particle conservation equation 
(and keeping only terms linear in the deviation away from 
equilibrium), we express the right side of Eq. (^) as 



(3.2) 



where Sv is the Eulerian velocity perturbation, and we 
assume that the perturbation has time dependence e~*"* 

in the comoving frame of the fluid. 

In order to analyze the left side of Eq. ( |3.l[ ), we exam- 
ine a fluid variable x that characterizes the microscopic 
process which produces bulk viscosity. For small depar- 
tures from equilibrium, the value of x in the physical fluid 
state relaxes toward its value in thermodynamic equilib- 
rium by 



dtx + V ■ Vx = — (a; — i)/r. 



(3.3) 



where r is defined as the relaxation time for this process. 
We are interested in nearly equilibrium fluid states in 
which the physical values of the state variable x (and 
hence the thermodynamic state x) oscillate about the 
background equilibrium, so that 

{dt + V ■ '^){x — xo) — ~iLj{x — xo), (3.4) 
(9t -I- w • V)(i — xo) = —i(l>{x — XQ). (3.5) 

In such a state it is straightforward to verify that 

X — X X Xq 



X — Xq 



lUJT 



1 — i6jT 



(3.6) 



Now consider how the fluid variable x changes as the 
particle number density of the state is varied slowly from 
one equilibrium state to another: 

X — xo = -^{h — no) = —An (3.7) 
an an 



(since by definition n = n). It follows then that the 
difference between p and po is given by 



P~Po 



dp 
dn 



dp 
dx 



dx 
dn 



An. 



(3.8) 



A similar argument (now using Eq. |3.6| to relate x — xq 
and X ~ Xq) gives the following expression for p — po'. 



P-Po 



dp 
dn 



1 





dx 


\dx^ 


'n dn_ 



An. (3.9) 



Combining Eqs. ( |3.8| ) and ( |3.9| ) gives us an expression for 
the difference in the pressure p — p that appears on the 
left side of Eq. (O: 



iujT dp\ dx ^ 
^ ^ 1 — i(2iT \ dx J ^ dn 



(3.10) 



Then equating this expression for the left side of Eq. ( |3.1| ) 
with the expression for the right side from Eq. (3.2), we 
find the desired formula for the bulk viscosity: 



c 



' dp 
1 — iiliT V dx 



dx 
dn 



(3.11) 



Finally it is convenient to re-expre ss th e thermody- 
namic derivatives that appear in Eq. (3.11) in terms of 
the more familiar 



n f dp 



loo = - [tt- 
p \on 



the "infinite" frequency adiabatic index, and 



n 

7o = - 
P 



dp 
dn 



dp\ dx 
dx J „ dn 



(3.12) 



(3.13) 



the "zero" frequency adiabatic index. In terms of these 
quantities, then the bulk viscosity may be written in the 
form: 



c 



p(7oo ~7o)t 
1 — Out 



(3.14) 



For a fluid composed of particles with non-relativistic 
energies this expression is equivalent to the conventional 
one [ p^ written in terms of the sound speed u, since 
P7 = pu^ for such fluids. However for a fluid contain- 
ing particles with relativistic energies, the conventional 
form is wrong and Eq. ( ^.14 ) is the appropriate form to 
use. Note that this new form of the bulk viscosity is 
needed to describe any fluid containing relativistic inter- 
nal particle energies, even if the bulk motion of the fluid 
itself has only non-relativistic velocities which arc well 
approximated by Newtonian hydrodynamics. We also 
see that for a fixed-frequency perturbation, the greatest 
bulk viscosity comes from processes with relaxation times 
T w l/o). The importance of this fact will become ap- 
parent as we examine how the relaxation time r varies 
inside a neutron star. 
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The standard approach jlj] to treating muhiple reac- 
tions is to repeat the preceding derivation for muhiple 
degrees of freedom Xi and relaxation times r^. However, 
this only works if (as in air) the Xi can be chosen to be 
independent. This is not possible in neutron-star mat- 
ter, since the degrees of freedom (e.g., concentrations of 
various baryons) are related to each other even out of 
thermodynamic equilibrium by constraints such as con- 
servation of baryon number. The reactions of interest to 
us here are the non-leptonic weak interactions 



n + n ^ p 



A. 



Given the microscopic reaction rates for these processes 
(which are calculated in Sec. IV), we can express all of 
the perturbed quantities in terms of a single one. Since 
all the hyperon reactions that contribute to bulk viscosity 
involve neutrons, we choose as our primary variable the 
number density of neutrons n„. 

Let Xn = Tin/nB be the fraction of baryons in a given 
fluid element that are neutrons. This variable changes 
only by internal reactions, not directly by changing the 
volume of the fluid element, and so we can write 



-{Xn- Xn)/T = -Tn/UB- (3.15) 



Here r„ is the production rate of neutrons per unit vol- 
ume, which is proportional to the overall chemical poten- 
tial imbalance 5^ = p — p (see below). (We normalize 
n„ and r„ to the baryon number density ub to remove 
the oscillating time dependence of the volume of the fluid 
element.) Assuming the reactions described in Sec. 
the relaxation time is then given by 



2r. 



Sfi UB^Xn 



(3.16) 



where Sxn = Xn — Xn- 

We obtain Spi/Sxn from the following constraints: 



= Sxn + SxA + Sxp + 6xy: , 
= Sxp — Sx^, 

= (3nSXn + f^A^XA + PpSXp 



(3.17) 
(3.18) 
/3sfeE, (3.19) 



where the number densities of baryon species nor- 

malized to the total number density of baryons and the 
f3i are defined below. The first constraint is conservation 
of baryon number, obeyed by all reactions. (We note 
that since Sub = nB — ub = that Sxi = Sui/uB-) 
The second constraint is related to conservation of elec- 
tric charge, but is stricter: We assume that all leptonic 
(Urea) reaction rates are much smaller than those which 
produce hyperon bulk viscosity and so protons are only 
produced in reactions that produce a E~. The third con- 
straint is that the non-leptonic reaction 



A 



P 



(3.20) 



proceeds much faster than the weak interactions which 
produce hyperon bulk viscosity since it is mediated by the 
strong nuclear interaction. Here we use the shorthand 



dpi 
9n, 



''3 / nk,k^j 

Pi — Oini + aAi — a. 



pi 



"Si 



(3.21) 
(3.22) 



Equilibrium with respect to reaction (B.20) ensures that 
both processes described in Sec. IV have the same chem- 
ical potential imbalance. 



5p = Spin - SpA = ^Spn - Spp - Sps- 



(3.23) 



It is straightforward then to express 6p in terms of the 
Sxi, and then to eliminate all but Sxn using the con- 
straints Eqs. ( |3.17D -( |3T9| ). The resuh, 



Sp 



iPn - PA){anp - OCAp + O-nY. - Q^As) 



ubSx^i 



-a An 



2/3a - /3p - /?s 
(2/3„ - I3p- l3T.){anA - oi aa) 
Wa -Pp-Pi: 



(3.24) 



then determines via Eq. ( 3.16| ) the relaxati on ti me that 



appears in the bulk viscosity formula Eq. (3.14). For a 
certain range of densities there are S~ hyperons present 
in Glendinning's equation of state, but no A. In that 
case the variable 5xa remains zero, and the constraint 
Eq. ( 3.19| ) is no longer enforced. In this case the chemical 
potential imbalance can still be expressed in terms of 5xn 
with the somewhat simpler result: 



25p 

ubSxt, 



4a„„ - 2{apn + a^n + cunp + a„s) 



asp + OLpY. + ass- 



(3.25) 



The equation of state of neutron-star matter is gen- 
erally written in a form which gives the thermody- 
namic variables as functions of the various particle 
species present, e.g. the pressure would be specified 
as p — p{ni). In order to evaluate the thermodynamic 
deriv ative s {dp/dn)^ and {dp/dx)n that are needed in 
Eq. (3.11), we note that n„ = nsXn for our choice of 



X. Thus the partial derivative needed in Eq. (3.11) is 
given by {dp/dx)n = nsdp/ drin. Similarly, if needed, 
{dp/dn)x = y^,- Xjdp / drij . The derivative dx/dn that 



also appears in Eq. (3.11) is determined by construct- 
ing a sequence of complete equilibrium models (e.g. by 
imposing all of the necessary /3-equilibrium constraints, 
etc.) for different total baryon number densities n^. This 
complete model of the equilibrium states will include the 
functions XnijiB), from which the derivative dxn/duB is 
easily computed. 

Formation of a superfluid (of a given particle species) 
is marked by the formation of Cooper pairs and col- 
lapse of the pairs into a Bose-Einstein condensate. It 
is the unpaired particles that we are concerned with, 
since they are the ones that can participate in the bulk- 
viscosity generating reactions. The free-particle states 
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within the pair-binding energy A of the Fermi surface 
are depleted. As a result all phase-space factors (and ef- 
fectively the reaction rates) are decreased by roughly a 
factor e~^/^'^ . The effect of superfluidity can be included 
in our ordinary-fluid bulk viscosity calculation then sim- 
ply by making the substitution 



fifth Dirac matrix is 7 = 17 7 7 7 , the me tric has neg 



(3.26) 



Thus when superfluidity is taken into account the equa- 
tion for the relaxation time, Eq. ( 3.16D becomes: 



(3.27) 



IV. MICROSCOPIC REACTION RATES 

Since the and A hyperons form at the lowest thresh- 
old densities, we are most interested in the nonleptonic 
reactions forming them from neutrons. Following Jones' 
earlier work [pi M we calculate rates for two reactions. 



n + n *^ p 



A, 



(4.1) 
(4.2) 



as tree-level Feynman diagrams involving the exchange 
oi aW boson. In his most recent paper Q, Jones treats 
the reaction 



■A, 



(4.3) 



which is generally the dominant channel for A produc- 
tion in laboratory experiments on hypernuclei. This pro- 
cess has no simple M^-exchange contribution. Several 
other processes contribute, but at present the rate cannot 
be well predicted from theorv Some o f th ese processes 
surely operate in reactions (4.1) and (4.2) and modify 
the rates, perhaps significantly. However, our simple cal- 
culations should provide a reasonable lower limit on the 
rates and thus an upper limit on the bulk viscosity. 



A. Single reactions 

We calculate reaction rates using the standard tech- 
niques of time-dependent perturbation theory in rela- 
tivistic quantum mechanics. We use the conventions of 
Griffiths pM: spinors are normalized to uu — 2m, the 



ative trace, and h — 1. For a single reaction (4T) or (4.2) 
between particles with 4-momenta pi , the differential re- 
action rate (number per unit volume per unit time) is 



dr = IX |2 (2^)4^(4) (pi +p2^P3-pi)Sl[ ■ 



^^l(27r)326,' 
(4.4) 

where is the spinor matrix element (squared and 

summed over spin states), boldface Pi are 3-momenta, 
and ei are particle energies. The statistical factor S, 
which compensates for overcounting momentum states of 
indistinguishable particles, is 1/2 for reaction (4.1) and 
1 for reaction (4.2). 

First, consider the A reaction ([l.2|), which we represent 
by a single tree-level diagram. Labeling the particles 1 
(neutron), 2 (ingoing proton), 3 (outgoing proton), 4 (A 
hyperon), we obtain the matrix element (for a single set 
of spin states) 



Ma = 



G 

-^-j= sin 26*0 ["(^3)7^ (1 + 5«p7^) u{Pi)u{P4) 
X 7m (1 + 5M7^) u{p2)] ■ (4.5) 



Here is the Fermi coupling constant and dc is the 
Cabibbo weak mixing angle. The quantities Qnpi 9nYii 
and are axial- vector couplings (normalized to the vec- 
tor coupling) of the weak interaction, whose deviation 
from —1 represents the partial nonconservation of the 
axial current. (These quantities are often written qa or 
Ga/Gv in the particle and nuclear physics literature. We 
add a label to keep track of which nuclcon-hyperon line 
is which.) We use the values Gp = 1.166 x lO^^MeV"^ 
and sin^c = 0.222 from the Particle Data Group [ pT| . 
The axial- vector couplings change with varying momen- 
tum transfer and density of the medium in a way that 
reflects the internal (strong-interaction) structure of the 
baryons and is therefore difficult to calculate. We use 
the values gnp = —1.27, — —0.72, and gnT, — 0.34 
measured in /?-decay of baryons at rest [|7| . There are 
theoretical reasons to believe that all axial-vector cou- 
plings tend to their asymptotically free values of — 1 in a 
dense medium We provide dissipation numbers in 

later sections both for laboratory values and for asymp- 
totic values of the couplings to give an estimate of the 
uncertainty in this calculation. 



To obtain the net reaction rate F we sum over all possible initial and final spin states. This is done in the 

standard way by tracing over outer products of spinors: 

\Mk\^ = AGp sin^ 200' {2m„TO^mA (l - gl-p) (l - gl^) ~ ™«mpP2 ■ P4 (l - glp) (l + 5pa) 
-mpiTLAPi • P3 (1 + glp) (1 - .9pa) + Pi • P2 P3 • P4 [(1 + glp) (l + .9pa) + 4g„p5pA] 

+Pl ■P4P2- P3 [(1 + glp) (1 + ffpA) - 45np5pA] } ' (4.6) 
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which in the low-momentum limit used by Jones Q reduces to 

\Ma\^ - 8G^p sin^ 20cmnmlmA (l + Sglpgl^) ■ (4-7) 

Note that this limit corresponds to Ga = OAQGp in the notation of Jones who obtains Ga = 1.29Gf- The 
factor of three discrepancy is within the uncertainties of modern nuclear-matter physics, as we see by e.g. taking the 
asymptotic values of the axial- vector couplings. 

Reaction (4.1) is treated similarly, with particle labels 1 and 2 (neutrons), 3 (proton), and 4 (E^ hypcron). Anti- 
symmetrizing with respect to the two indistinguishable neutrons, the matrix element is 

G F 

Ms = ^ sin 26lcu(p3)7^ (l + 5«p7^) [w(pi)m(P4)7m (l + 5nS7^) "(^2) - u{p2)u{p4ht. (l + ffnS7^) ■ (4-8) 

The squared sum over spins is given by 

\Msf = 4:Gp sin^ 200 {em^TOpms (l - g^p) (l - gl^) - ■mpmj:pi ■ P2 {l ~ glp) (l - Sns) 
-2m„mspi • Pa (1 + glp) (l - gl^) - 2m„mpPi ■ P4 {l - glp) (l + gl^) 
-2rn„msP2 • Pa (l + glp) (l - gls) - 2rn„"ipP2 ■Pi{l~ glp) (l + .9^s) 

-mips ■ p4 [(1 -I- glp) (1 + ,9,^5.) - 4.g„p,g„s] + 4pi • P2P3 • P4 [(l + (l + ff^s) + 4.g„p5„s] 

+Pi ■P3P2 ■P4 [(l+.g^p) (l+.gfjE) -4.9„pgns] +Pi ■P4P2 -Pa [{^ + glp) (1 + .9^e) -4g„p5„s] }, (4.9) 



which in the low-momentum limit reduces to 



where 



lAlsP = SG% sin^ 20cmlmpms (1 -f 3.9„p5„s) • 

(4.10) 

This expression agrees with Eq. (7) of Jones' old pa- 
per Q, allowing for our different placement of the sta- 
tistical overcounting factor and different normalization 
of the spinors (he uses uu = 1). However, we find that 
for this reaction the low-momentum limit is a very poor 
approximation to the full result: the collision integral 
computed in the next subsection c an be more than an 
order of magnitude higher than Eq. ( 4.10 ) would suggest. 
This effectively brings down the coefficient of bulk viscos- 
ity C by a factor of about 5 at densities 5-6x10^'' g/cm"^, 
where there are no A hyperons MM . 



B. Collision integrals 



Reactions (4.1) and ( [4.2| ) can be regarded as scattering 
processes wherein the scattered particles change identity. 
Therefore we can use existing results on collision integrals 
in the literature on superfluidity (e.g. |po|| ) with some 
slight modific ation s. We now integrate the differential 
reaction rate (4.4) over momentum space to obtain the 
total rate 



S 



t[^\M\'S('Hpi 



P2 P3 -P4) 

(4.11) 

(From here on we use italic pj to denote the absolute value 
of the 3-momenta, pi = y/Pi ■ Pi ) The Pauli blocking 
factor, 



40967r8 

y.F{ei)5{ti + e2 - - ^4) 



/, = l/{l + exp[(e,~A^,)AT]}, 



(4.13) 



accounts for the degeneracy of the reactant particles 
and restricts the available phase space to those particles 
within roughly kT of their Fermi energies. 

In the case where all particles are degenerate (which is 
true except for a very small region just above the thresh- 
old density for each hyperon species), the collision in- 
tegral separates into angular and energetic parts. The 
energetic part can be written in the limit 5^ <C kT as 



j ]Jde,F(ei)^(ei + £2 - £3 - £4) 

i 

y'^dy 



{kTySfi 



+ 00 



{ev - 1) (1 - e-y) ■ 



(4.14) 



where the latter integral has the value 27r^/3. 

We address the angular integral with the aid of Fig. |^. 
With no preferred direction (e.g., strong magnetic field), 
the pi volume element is Anpfdpi and the p2 volume ele- 
ment is p\dp2 sin OdO d4>, where 6 and 4> are the polar and 
azimuthal angles of P2 with respect to pi . Conservation 
of momentum demands that pi -\- P2 = Ps + P4, deter- 
mining a common axis. Let a and a' be the angles of 
Pi and P3 with that axis. The volume element of pa is 
then p3 sina'd0', where (f)' is the angle between pi x p2 
and p3 X p4, multiplied by the area element in the plane 
containing p^ and P4. To separate the energetic and an- 
gular integrals, it is convenient to write this area element 
as dp^dp4/ sin 6'', where 0' is the angle between p3 and P4. 
Before leaving Fig. 0, note the following useful identities: 



F{e,) = /i/2(l - /3) (1 - /4) - (1 - /l) (1 - /2)/3/4, (4.12) 



2p\p2 cos Q 



2 
P3 



pl + 2p3P4COse', (4.15) 
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FIG. 7: Definition of angles used in the collision integral. 
The plane containing pa and p4 can be rotated out of the 
page around the long-dashed common axis by an angle 0' for 
a given pi and P2 • 



Pisina = P2sin(0 — a), (4.16) 
Pasina' — p4 sin(0' — a'). (4-17) 

We use the differential of the first identity (with -pi held 
constant) to write our final result for the volume element 
as 



47rp3sina'd6''d0'd0 J|dej. (4.18) 



Although (/) and 0' are free to range from to 27r, 
the limits of integration over Q' depend on the relations 
between the momenta pi, which are constrained to be 
close to the Fermi momenta. (Therefore from now on we 



use Pi to refer to the Fermi momenta.) We integrate d' 
over the full range to tt, which is allowed by momentum 
conservation for 

Pn>Pp>iPA,P^), Pn-Pp<Pp-PA- (4.19) 

In our equation of state these Fermi momentum criteria 
are satisfied for p < 1.0 x lO^^g/cm'^ (see Fig. |l). It turns 
out that almost all of the dissipation takes place below 
this density, and much above this density the baryons are 
too closely packed to maintain their separate identities 
anyway. Therefore for our purposes it is sufficient to 
treat only the case ( 4.19| ). 

Now consider the case that does not depend on 

any of the angles, such as in the limit Pi/rrii ^ 1. We can 
integrate trivially over all the angles, with the exception 
of 



/ p3 sin a' d9' — 2p4, 
Jo 



(4.20) 



where p4 is the smallest (i.e. hyperon) Fermi momentum 
involved. This motivates us to use Eq. ( 4.14 ) to write 



r 



S 



-{\M\^)p4kT)^Sfi, 



, 1927r3 „ 
where (|A^| ) is the angle-averaged value of 

write |A^p in terms of the integrals 



(4.21) 
We 





P2> 


= i{~3pl~3pl + 3pl+pl), 


(4.22) 




Pa) 


= i{-3pl + 3pl-3pl+pl), 


(4.23) 


(Pl 


P4) 


= (P2 • Pi) = (P3 • Pi) = -ipl, 


(4.24) 




P3> 


= U^pI-^pI-SpI+pI) , 


(4.25) 


Pi • P2 P3 


Pi) 


= ^pK^pI + ^pI + ^pI-pI), 


(4.26) 



(Pi • P3 P2 • Pi) = (Pi ■PiP2-P3) 



Pi 



60 (pI 



pI) 



5pt - IOpIpI + 5pI + 3 {5pt - GpIpI + pt) 



(4.27) 



(Note that some dot products, e.g. pi p^, depend on (p' . The average over (p' is then nontrivial but is easily taken by 
symmetry about the common axis.) 
The results for the angle averages are 



G% sin^ 29c 
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{120 (1 - glp) (1 - m„m^mA - 20 (l - gl^) (l + g^^) mnUip (SeptA - pI) 



-10 (1 + glp) (1 - 5^a) mpUiA (6e„ep - 3pl + pi) + 2 [(l + gl^) (l + .g^^) + ig^pgpA] 

X [5epeA (6e„ep 4- 3pl - pi) + pi (l0e„ep + 5pl + lOp^ - pi)] + [(l + g^p) (l + g^^) - 4g„p5pA] 

X [lOe„eA {6ml + 3pl + pi) + pi {-20el + Ibpl ~ 3pl + h{pl - p^Vipl - pI))] } . (4.28) 



[Note that the denominator in the last term does not diverge while criterion (4.19) holds.] Similarly, for reaction (4.1), 



{\M, 



= Y^G% sin^ 2ec {180 (l - g^p) (l - gl^) mlmpnis - 40 (l - glp) (l + gl^) m„mp (3e„es - p|) 
-20 (l + glp) (l - glj:) m„ms (6e„ep - 3p^ -I- p|) - 5 (l - glp) (l - glj^) mpniT. (6e^ -I- 6pl - 3pl 



pD 



10 



(1 



.9ns) 



+4[(1 

+ 10 [(1 + glp) (1 + .g^j.) - 4g„pg„E] [-m^ (Stpes + p|) + e„ (6e„epes - 2epp| - 3e^pl + esp|)] } ■ (4.29) 



These angle averages are inserted into Eq. (4.21) to yield 



the net reaction rate per unit volume as a function of 
Fermi momenta. 



V. RELAXATION TIMES 

The relaxation tim escale s fo r the non-leptonic weak 
interactions, Eqs. (4.1) and (4.2), can now be determined 
by combining the microscopic collision rates determined 
in Sec. IV with the thermodynamic quantities evaluated 
in Sees. II and III. At densities above the threshold for 
the production of S^, but below the A threshold, the 



final expression for the relaxation timescale r is, 



{kT)^p^{\Mj 



1 



r 1927r3 e'^^/'^^ nsSxri' 



(5.1) 



Here the c ollision cross section |A4sP is evaluated from 
Eq. ( 4.29 ). The Fermi momenta that appear in this 
expression are obtained from the complete description 
of the equilibrium thermodynamic state as described in 
Sec. II. This description includes the values of the parti- 
cle densities (and hence Fermi momenta) of each species 
as a function of the total baryon density. The thermo- 
dynamic quantity dfi/nB5xn is given by Eq. (3.25) in 



terms of chemical potentials and their derivatives. Once 
the density increases to the point that both E~ and A 
hyperons are present in the equilibrium state, then the 
expression for the relaxation timescale becomes, 



1927r3 



, pa{\Ma 



2\1 



oAA/fcT 



(5.2) 



Here the collision cross sections are given by Eg . ( 4.2S| ) 
and ( 4.29| ), and Sii/nsSxri is given by Eq. ( 3.24 ) in this 
case. 

We have evaluated this relaxation timescale for 
neutron-star matter using the equation of state described 
in Sec. II. Figure ^ illustrates the density dependence of 
this timescale for a range of temperatures. In the case 
of Fig. H we have assumed that the S~ superfluid gap 
function is given by As — Aa, while in Fig. ^we assume 
As — 10 Aa. The only significant difference between 
these two cases comes about in the density range where 
there exist but not A. In that range the timescale 
is significantly increased in the As — IOAa case by the 
stronger superfluid effects. All of the curves in these two 
figures were evaluated using the "standard" /3-decay val- 
ues of the axial-vector coupling constants: Qnp — 

-1.27, 

5pA = —0.72, and 5„s = 0.34 |17 and the effective 
masses of all of the baryons when evaluating the scatter- 
ing cross sections. There are theoretical arguments that 
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FIG. 8: Density dependence (in units of 10^^ g/cm'^) of the 
relaxation timescale r (in units of s) for a range of temper- 
atures. These curves were constructed using the assumption 
Ay^- = Aa for the E~ superfluid gap. 
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FIG. 9: Density dependence (in units of lO'^^ g/cm"^) of the 
relaxation timescale r (in units of s) for a range of temper- 
atures. These curves were constructed using the assumption 
Aj^- = IOAa for the E~ superfluid gap. 



suggest the coupling constants should approach the val- 
ues, gnp — 5pA = ffnS = — 1 in a dense medium jl^. We 
illustrate the impact this might have on these timescales 
in Fig. |l^ (dotted curve). We also illustrate in this fig- 
ure (dashed curve) the effect of using the bare masses of 
the various baryons when computing the scattering cross 
sections. We see that the overall effect of these changes 
is to make the timescales shorter (by up to an order of 
magnitude). This tends to decrease the bulk viscosity by 
a similar factor, until the temperature drops below the 
superfiuid critical values. 

Finally, we are in a position now to evaluate the bulk 
viscosity itself. The real part of the bulk viscosity, the 
part that is responsible for damping the modes of neutron 
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32 



FIG. 



10: Density dependence (in units of 10^^ g/cm^) of the 



relaxation timescale r (in units of s) for T — 10 K. Solid 
curve uses effective masses and the /3-equilibrium values of the 
axial-vector couplings. Other cuves explore various alternate 
microphysics assumptions: dashed curve uses bare masses, 
dotted curve uses gs ~ — 1- 



stars, is given by the expression 

Piloo - 1o)t 



Re C 



1 + {uJTf 



(5.3) 



Figure |Tl| illustrates the density and temparature depen- 
dence of C,. (Here we assume that the frequency corre- 
sponds to the m — 2 r-mode frequency of a maximally 



rotating neutron star: uj 



r^^max-) These figures illus- 



trate the complicated temperature dependence of the vi- 
socity due to superfluid effects. For temperatures slightly 
below the superfluid critical temperature the values of 
the bulk viscosity are increased over their normal val- 
ues. This is due to an increase in the timescale r which 
moves it closer to being in resonance with the pulsation 
period of the mode. Once the temperature falls well be- 
low the superfluid critical temparature however, we see 
that the timescale r becomes even longer than the pul- 
sation period and so the viscosity becomes smaller again 
in this case. We note that even for very low tempera- 
tures there exists a small range of densities, just above 
the hyperon threshold densities, where the bulk viscos- 
ity remains rather large. This is due to the momentum 
dependence of the superfluid gap, Eq. (p.?!). The gap A 
goes to zero as the Fermi momentum of the particle goes 
to zero. Thus just above the threshold density the su- 
perfluid gap vanishes (for any finite temperatue) so the 
material in this region will retain the normal-fluid value 
of the bulk viscosity. 

Our value of C is generally much larger than that ob- 
tained recently by Jones H: at a total density p — 

7 X lO^^'g/cm'^ and temperature T — 10^°K, our ( is 
larger than Jones' by a factor of 400. Roughly a factor of 

8 is due to the r elat ively weak coupling we calculate for 
reaction in Eq. (4.2). (At this density the E~ hyperons 



account for only about 10% of the bulk viscosity and can 
be neglected.) We note that using the asymptotic values 
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FIG. 11: Density dependence (in units of 10^^ g/cm^) of the 
hyperon bulk viscosity (in units of g/cm s) for a range of 
temperatures. These curves were constructed using the as- 
sumption Ag- ~ Aa for the E~ superfluid gap. 




FIG. 12: Density dependence (in units of 10^^ g/cm^) of the 
hyperon bulk viscosity (in units of g/cm s) for a range of 
temperatures. These curves were constructed using the as- 
sumption Aj;- = IOAa for the E~ superfluid gap. 



of the weak axial-vector couplings would erase much of 
this factor of 8 difference, and thus it is indicative of the 
size of the uncertainties in ^ due to our poor understand- 
ing of nuclear-matter physics. The remaining factor of 
50 is thermodynamic in origin. Jones evaluates various 
partial derivatives of pressure and chemical potentials, 
e.g. appearing in Eq. (3.24) of our paper and Eq. (42) 
of Ref. §, using the values for a gas of noninteracting 
fermions. We include all the mesonic interaction terms, 
whose effect is to increase significantly these thermody- 
namic derivatives. Since the details of the neutron-star 
equation of state are uncertain, the precise values of the 
derivatives are correspondingly uncertain. However, we 
think it unlikely that the true physical mesonic terms will 
cancel precisely enough to bring the derivatives down to 
their noninteracting values. In summary, we think that 
the true value of C is within an order of magnitude of the 
value we compute here. 
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VI. r-MODE DAMPING TIMES 

The basic formahsm for evaluating the effects of bulk 
viscosity on the stability of the r-modes is well known ^, 
pT| . The time derivative of the co- rotating frame energy 
E due to the effects of bulk viscosity is 



dE 
'dt 



Re C |V 



(6.1) 



This rotating frame energy E is (to lowest order in the 
angular velocity of the star) given by the integral 



E=- I p\Sv\'d^x. 



(6.2) 



The timescale ''"^(/i) on which hyperon bulk viscosity 
damps the mode then is 



1 _ 1 dE 
TB(h) 2E dt ■ 



(6.3) 



Here we have normalized T5(h) so that 1/T5(h) is the 
hyperon bulk viscosity contribution to the imaginary part 
of the frequency of the mode. 

For the case of the r-modes (in slowly rotating stars) 
the i nteg rals that determine E and dE/dt in Eqs. (6.1) 
and (5.2) can be reduced to simple one-dimensional inte- 
grals. For the case of E this reduction is well known M: 



E = -a^n^R-^ 
2 



pr dr. 



(6.4) 



Here a represents the dimensionless amplitude of the r- 
mode, and Q and R are the angular velocity and radius 
of the star respectively. The reduction of dE/dt to a 
one-dimensional integral is not so straightforward. In 
general the expansion of the mode V ■ 5v is a compli- 
cated function of radius and angle. To lowest order in 
slowly rotating stars the bulk viscosity C, will de pend only 
on radius. Thus we may always convert Eq. (6.1) to a 
one-dimensional integral by defining the angle averaged 
expansion squared (|V • 5v\'^): 



dE 
'dt 



= -47r 



Re C (|V • Sv\^ 



r^dr. 



(6.5) 



While the angle-averaged expansion is in general a com- 
plicated function, for the case of the r-modes it is rather 
simple. This function has only been determined numeri- 
cally pH], however the simple analytical expression, 



(6.6) 



is an excellent fit to those numerical solutions. 

Once the structure of the density function p{r) in a 
stellar model is determined, it is straightforward to eval- 
uate the integrals in Eqs. (|6.4D and (6.5) using Eq. (|6.6|) 



to determine the bulk viscosity damping time ts(/i) . The 
bulk viscosity of interest to us here is very sensitive to the 
density of hyperons in the stellar core. Thus we use the 
relativistic stellar structure equations to evaluate p{r). 
This ensures that the size and stucture of the hyperon 
containing core are sufficiently accurate for our purposes. 
These functions p{r) are illustrated for a range of stellar 
masses in Fig. ^ based on the equation of state discussed 
in Sec. II. Given these (numerical) expressions for p(r) 
it is straightforward then to use the expressions for the 
hyperon bulk viscosity C derived in Sec. ^ to obtain C(r) 
for any given neutron star temperature. Together p{r) 
and C(r) then determine TB(h) through Eqs. (|6.4|), ( 6^ ) 
and ( |6.3| ). While it is straightforward to evaluate these 
timescales, the result is a rather complicated function of 
the temperature, angular velocity, and mass of the stellar 
model and so we do not attempt to illustrate it directly. 

The most important application of the hypcron bulk 
viscosity timescale ''■^(/i) is the analysis of the role this 
type of dissipation plays in the gravitational radiation 
driven instability in the ?'-modes. Gravitational radia- 
tion contributes a term to the evolution of the energy 
dE/dt that is positive. As is well known by now, gravi- 
tational radiation tends to drive the r-modes unstable in 
all rotatingstars ^ . As has been discussed in detail 
elsewhere |£[ ^ the imaginary part of the frequency of 
the r-mode may be written as: 



1 



1 



1 



1 



(6.7) 



Here tq]^ represents the timescale for gravitational radi- 
ation to effect the r-mode, TB{h) is the hyperon bulk vis- 
cosity timescale discussed here, and Tb{u) is the modified 
Urea bulk viscosity. Detailed expressions for evaluating 
these other terms are discussed elsewhere and will not 
be repeated here. Suffice it to say that each is a func- 
tion of the temperature, angular velocity and mass of the 
neutron star. Since l/r^ is the imaginary part of the fre- 
quency of the r-mode, the mode is stable when > 
and unstable when < 0. For a star of given tempera- 
ture and mass, the critical angular velocity D.^ is defined 
to be the angular velocity where I/t^ = 0. Stars rotating 
more rapidly than f2c are unstable while those rotating 
more slowly are stable. 

We have evaluated the critical angular velocities fie 
numerically using the new hypcron bulk viscosities de- 
rived in Sec. ^. Figure 13| illustrates the temperature 
dependence of the critica angular velocities for a range 
of neutron-star masses. The more massive neutron stars 
have larger hyperon cores which suppress the r-mode in- 
stability more effectively. The curves in Fig. |l^ assume 
that the superfluid gap function is given by As = Aa, 
and that the axial vector coupling coefficients have their 
/3-decay values. In Fig. |l^ we compare the critical angu- 
lar velocity curves for 1.4M0 stellar models using either 
the As = Aa (solid curve) or the As = 10 Aa (dash-dot 
curve) assumption about the E~ superfluid gap. The 
larger value of As allows superfluidity to make the bulk 
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FIG. 13; Critical angular velocities for neutron stars as a 
function of hyperon core temperature. Each curve represents 
a neutron star of fixed mass, ranging from 1.2 Mq to the 
maximum mass for this equation of state, 1.79 Mq. These 
curves assume As = Aa and use the /3-decay values of the 
weak axial vector coupling coefficients. 



FIG. 15; Critical angular velocities for IAMq neutron stars as 
a function of hyperon core temperature. Solid curve includes 
the effects of hyperon bulk viscosity, hyperon channel direct 
Urea bulk viscosity, and modified Urea bulk viscosity. Dashed 
curve leaves out the effects of the direct Urea hyperon bulk 
viscosity. 




FIG. 14; Critical angular velocities for 1.4 Mq neutron stars 
as a function of hyperon core temperature. The solid curve 
assumes that Ae = Aa while the dot-dash curve assumes 
Ae = 10 Aa. Both curves use the /3-decay values of the weak 
axial vector coupling coefficients. Dotted curve uses gs ~ —1, 
while dashed curve uses bare masses. 



viscosity larger over a wider range of temperatures, and 
hence the r-mode instabiUty is less effective. Also il- 
lustrated in Fig. |l^ are the effects of changing various 
microphysics assumptions. The dotted curve shows the 
effect of changing the values of the axial vector coupling 
constants from their /3-decay values to the asymptotic 
value -1. And the dashed curve shows the effect of us- 
ing bare masses rather than the effective masses in the 
scattering cross sections. These alternative assumptions 
make the bulk viscosity less effective and the r-mode in- 
stability operates over a wider range of angular velocities 
in these stars. However, neither of these effects is as large 
as that resulting from a change in the superfluid gap. 
The hyperons' primary contribution to the bulk vis- 



cosity of neutron-star matter is through the mechanism 
discussed above. However, as pointed out by Jones, the 
presence of hyperons in the core of a neutron star also 
make it possible for alternate forms of the direct Urea 
interaction to take place and these too contribute to the 
bulk viscosity of the material. Jones showed that the 
contributions to the bulk viscosity from this process are 
given by 



Rc C = 



4.9 X IO^OTi 



1 



2.0 X lO-'^oj'^T^Q^ ' 



(6.8) 



in cgs units for typical values of neutron star matter, 
where Tio is the temperature measured in units of 10^'^K. 
Using this expression in the A containing core of the neu- 
tron star, we have evaluated the effects of this hyperon 
channel direct Urea bulk viscosity on the stability of the 
r-modes. These results are illustrated in Fig. |l^. The 
solid curve in Fig. |l^ includes the effects of the hyperon 
bulk viscosity discussed above, t he h yperon channel di- 
rect Urea bulk viscosity of Eq. (6.8), and the ordinary 



modified Urea bulk viscosity. For comparison the dashed 
curve leaves out the effects of the hyperon channel direct 
Urea bulk viscosity. We see that this direct Urea bulk 
viscosity has only a small effect on the stability of the 
r-modes for temperatures around lO^^K. 



VII. DISCUSSION 

We have analyzed here the effects of the bulk viscos- 
ity due to hyperons on the stability of the r-modes in 
rotating neutron stars. Hyperons exist only in the high 
density core of a neutron star where the influence of the 
r-mode is quite small. Thus to evaluate accurately and 
reliably the importance of this effect, it was necessary 
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to compute detailed and accurate models of the com- 
position and structure of the neutron star core, and to 
have an accurate model of the structure of the r-mode 
in this region. We use Glendenning's jl^, |l^ relativistic 
mean-field equation of state to evaluate the composition 
of the nuclear matter in the stellar core, and solve the 
relativistic Oppcnheimer-Volkoff equations to determine 
the stellar structure. 

Our evaluation of the hyperon bulk viscosity improves 
on previous work in se veral ways. First we generalize in 
Eqs. (3.11) and (3.14) the standard expression for the 
bulk viscosity coefficient so that it applies to relativistic 
fluids such as neutron st ar m atter. Seco nd we generalize 
the expressions in Eqs. ( 3.24 ) and ( 3.25 ) for the thermo- 
dynamic quantities that relate the microscopic reaction 
rates to the relaxation time that appears in the expres- 
sion for the bulk viscosity by including fully i ntera cting 
nuclear matter. And third we obtain in Eqs. (4.28) and 
(4.29) the fully relativistic expressions for the relevant hy- 
peron scattering cross sections needed to evaluate the mi- 
croscopic reaction rates. While our expressions for these 
cross sections reduce to the published low- mom entum re- 
sults, we find that the difference of Eq. ( 4.29 ) from the 
low-momentum limit can be an order of magnitude or 
more and reduces the coefficient of bulk viscosity some- 
what at low densities. Finally we evaluate the effects of 
this hyperon bulk viscosity on the r-modes using an ac- 
curate model for the structure of the r-mode in the core 
of a neutron star. 

Our results show that the hyperon bulk viscosity does 
not substantially suppress the gravitational radiation in- 
stability of the r-modes until the temperature of the core 
of the neutron star drops below a few times lO^K. (This 
is in spite of the fact that our coefficient of bulk viscosity 
is actually higher than that of Jones ||, ^ . The expan- 
sion of the fluid in the core of the star as given in Eq. 6.6 
is smaller than he estimated.) Below lO^K the r-mode 
instability is strongly suppressed in all of our models over 
the range of S~ superfluid gap functions and the range 
of axial vector coupling constants that we studied. If 
the core of the neutron star cools according to the stan- 
dard modifled Urea process |^ , then it would remain hot 
enough for the r-mode instability to act for about a day. 
This is enough time for the r-mode to grow and radiate 
away through gravitational waves a substantial fraction 
of the rotational kinetic energy of a rapidly rotating neu- 
tron star . However if the core of the neutron star cools 
substantially faster than this, then it may not be possible 
for the r-mode to grow rapidly enough to effect the star 
in a substantial way before the hyperon bulk viscosity 
stabilizes it. Cooling by the direct Urea process is signif- 
icantly faster that the modified Urea process: cooling the 
core of a neutron star to a few times lO^K within about a 
second |8| . Cooling by the direct Urea process will oc- 
cur in neutron-star matter whenever the proton/baryon 
ratio is larger than about 0.15. Since proton fractions in 
excess of this are now generally expected in neutron star 
matter, cooling by the direct Urea process seems likely 



at least until the temperature of the core falls below the 
superfiuid transition temperature for neutrons or protons 
at about lO^K. Thus it appears likely that the r-mode 
instability is effectively suppressed by rapid cooling of 
the neutron star core and the non-leptonic hyperon bulk 
viscosity. 

Once a neutron star cools below the transition temper- 
ature for the formation of neutron and proton superfiuids, 
the relaxation timescale for the hyperon interactions will 
increase exponenti ally compa red to the expressions de- 
rived here in Eqs. (5.1) and (5.2). This sharply reduces 
via Eq. (5.3) the bulk viscosity from this process at suf- 
ficiently low temperatures. Further detailed calculations 
would be needed to determine whether the hyperon bulk 
viscosity has a significant influence on the r-mode in- 
stability at temperatures of a few times 10* K, which are 
expected to exist in the cores of neutron stars in low-mass 
x-ray binaries. We did not carry out these calculations in 
part because solid crust-related shear |2^ and mag- 
netic field 1^ effects are quite effective in suppressing 
the instability at these low temperatures. 

How robust is the conclusion that the r-mode instabil- 
ity is effectively suppressed? Clearly the details of the 
nuclear physics involving hyperons in neutron star mat- 
ter are not well understood at this time. However, our 
conclusion applies to the entire expected range of the 
most poorly known properties of this material: the su- 
perfluid gap function and the axial vector coupling 
coefficients. In order to escape this conclusion, it would 
be necessary for neutron star matter to have very few 
hyperons present at the densities which exist in the cores 
of real neutron stars. This would require the equation 
of state to be substantially different from the one stud- 
ied here, or the masses of neutron stars to be signifi- 
cantly smaller than 1.4M0. Rapid rotation also lowers 
the central density and consequently the size of the hy- 
peron core in a neutron star. The central density of a 
maximally rotating 1.4M0 neutron star is about 73% of 
its non-rotating value (for the equation of state studied 
here) p7|| . This reduction almost eliminates the hyperon 
core for this extreme angular velocity, but over almost all 
of the range of angular velocities, 1.4M0 stars have sub- 
stantial hyperon cores. Finally if the dissipation in the 
core were sufficiently large it might be possible for the 
r-mode eigenfunction to be clamped to zero in the core 
by the dissipative processes while remaining finite and 
unstable in the outer parts of the star. The discussion of 
this possibility in the appendix shows that the hyperon 
bulk viscosity is not strong enough to clamp the r-mode 
in this way. 
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locally in one oscillation period is comparable to the lo- 
cal energy density of the mode. Thus we define the local 
quality factor of the mode: 



uie 
2ttI' 



(A4) 



APPENDIX A: MODE CLAMPING 



If g ^ 1 the mode will be clamped. For the r-modes we 
find that 



Bulk viscosity damps a mode by dissipating energy ac- 
cording to the expression 



dE 



E 

2^' 



(Al) 



For the case of the r-modes in slowly rotating stars, we 
may express the energy, and its time derivative as simple 
radial integrals: 



E : 



edr. 



dE 
'dt 
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(A3) 



where e and e are the angle averaged energy and en- 
ergy diss ipatio n ra te densities respectively as given in 
Eqs. Q and (fj). 

The mode will be completely suppressed (clamped) lo- 
cally if the amount of energy removed from the mode 
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For the IAMq neutron star model considered here 
^^max ~ 4700 rad/s, p > 5 x IQi^ g/cm^ and r < 6 
km in the region where hyperons occur, and i? w 14 km. 
Thus 



c 



8 X 1031 n 



(A6) 



for this case. Figures pi] and [T^ show that the bulk viscos- 
ity never exceeds about 10'^^ in the temperature range of 
interest to us. Thus, we conclude that the hyperon bulk 
viscosity represents a small perturbation on the basic hy- 
drodynamic forces of the r-modes. The condition g < 1 is 
violated only for relatively slowly rotating stars. In the 
domain where the gravitational radiation instability is 
most likely to be important, the dissipation by hyperons 
represents a small perturbation on the basic hydrody- 
namic forces, thus the r-modes will not be clamped. 
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